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Abstract 

We study the classical theory of a non-Abelian gauge field (gauge group SU{2)) coupled 
to a massive dilaton, massive axion and Einstein gravity. The theory is inspired by the 
bosonic part of the low-energy heterotic string action for a general Yang-Mills field, which 
we consider to leading order after compactification to (3 + 1) dimensions. We impose the 
condition that spacetime be static and spherically symmetric, and we introduce masses 
via a dilaton-axion potential associated with supersymmetry (SUSY)-breaking by gaugino 
condensation in the hidden sector. In the course of describing the possible non-Abelian 
solutions of the simplified theory, we consider in detail two candidates: a massive dilaton 
coupled to a purely magnetic Yang-Mills field, and a massive axion field coupled to a non- 
Abelian dyonic configuration, in which the electric and magnetic fields decay too rapidly 
to correspond to any global gauge charge. We discuss the feasibility of solutions with and 
without a nontrivial dilaton for the latter case, and present numerical regular and black hole 
solutions for the former. 



1. Introduction 



Following the early investigation of a variety of field theories coupled to Einstein gravity 
[1-4], it was widely believed that only the charges carried by massless gauge fields could 
characterize the exterior of a black hole. The notion that mass, angular momentum, and 
"electric" and "magnetic" charges are the only distinguishing features outside the horizon 
became known as the no-hair conjecture. 

In light of the no-hair results and several no-go theorems for classical glueball solutions 
with [5] and without [6] gravity, the recent discovery of both black hole [7-9] and smooth [10] 
solutions of SU{2) gauge theory coupled to Einstein gravity came as quite a surprise. The 
fields in such solutions decay sufficiently quickly that no global gauge charges are present, 
and hence no imprint at spatial infinity is required for the existence of nontrivial gauge field 
structure. It was later shown that static solutions with global electric or magnetic gauge 
charges can only occur in the embedded Abelian sector of this theory, and that non-Abelian 
dyons and dyonic black holes are prohibited [11]. Further analysis has also established a 
sphaleron interpretation of some smooth solutions which bridge topologically distinct Yang- 
Mills vacua [12-13], and the inherent instability of such saddle-point field configurations 
may help explain the generic instability of all Einstein- Yang- Mills (EYM) solutions against 
collapse into Schwarzschild black holes [14-17]. Despite their lack of stability, these non- 
Abelian solutions still present a challenge to the no-hair results, which are not based on the 
issue of stability, but rather follow from the careful analysis of several theories which are 
fundamentally different in character from EYM theory. 

In fact, the advent of these solutions has helped inspire a re-thinking of the no-hair 
conjecture, as well as a wealth of other solutions incorporating non-Abelian structure. In 
[18], a distinction is drawn between primary hair, such as the structure arising from the 
familiar continuous gauge charges, and secondary hair [19], which exists solely as a result of 
primary hair sources and hence docs not constitute a fundamentally new characteristic 
This distinction is well illustrated by the two recent approaches to Einstein- Yang-Mills - 
Higgs (EYMH) theory: for the "black holes inside magnetic monopoles" of [20-22], the 't 
Hooft-Polyakov monopole charge supports a core of secondary (triplet) Higgs hair outside 

#1 Wc do not discuss quantum hair in this paper, although this distinction applies equally well to the 
quantum and classical cases. 
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the horizon, while for the case of SU{2) coupled to a Higgs doublet (the massive vector 
theory of the standard model less hypercharge) [23], gauge and Higgs hair exist near the 
horizon without global gauge or topological charges. Thus, a distant observer in the latter 
case would not be able to distinguish such an object from a Schwarzschild black hole of the 
same mass, which motivates an alternate definition of primary hair: when the properties of 
a black hole are no longer completely determined within a given theory by the mass, angular 
momentum and continuous gauge charges, the additional parameters required to describe 
the black hole expand the space of states and give rise to primary hair [18]. 

Though the latter case and the original EYM black holes are examples of such primary 
hair, neither seem to share the stability properties of the well-known primary hair solutions. 
The black hole solutions of the spontaneously broken gauge theory appear to be unstable 
because of their similarity to the sphaleron-hke EYM solutions of [10, 7-9], and their inter- 
pretation as gravitating generalizations of the familiar SU (2) sphalerons: the weak-gravity 
hmit of one class of solutions is equivalent to the YMH configuration of [24-25]. On the 
other hand, when we ignore Hawking radiation, the secondary hair solutions of [20-22] are 
stable for the same reasons flat-space monopoles are stable. Thus, the physically important 
condition of stability appears to be more closely tied to the stability properties of correspond- 
ing flat-space solitons (when such solutions exist) than to the classiflcation of structure as 
primary or secondary hair. Other ilhistrations of this point have been found (e.g. the lin- 
early stable [26] black hole solutions to Einstein-Skyrme theory [27-28]), and a systematic 
approach to the existence of black hole solutions with solitonic flat-space counterparts has 
been formulated [29] , but a systematic treatment of the stabihty of such solutions is currently 
lacking. The search continues for stable black hole solutions to physically relevant theories, 
even several years after the flrst challenge to the no-hair conjecture opened the door to rich, 
new structure in the exterior of black holes. 

On a separate front, some recent progress in black hole physics has stemmed from the 
generic modiflcations to gravity mandated by string theory, a promising candidate for a 
consistent theory of quantum gravity which also provides predictions that challenge general 
relativity well below Planck scale curvatures. In particular, the presence at low energies 
of the dilaton and axion, two scalars with unusual couplings which appear in the same 
supersymmetric multiplet as the graviton, has precipitated a host of new black hole solutions 
with secondary hair and interesting properties. In charged dilaton black holes [30-35], the 
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Maxwell field acts as a source for dilaton hair, which leads to modifications of causal structure 
that help shed light on several puzzles peculiar to the Reissner-Nordstrom spacetime, as well 
as some mysteries of the later stages of Hawking evaporation. Because the axion couples 
to FF ~ E • B, black holes with both electric and magnetic charge can support axion hair 
[36,37]. Another axion coupling is of the Lorentz Chern-Simons form, so that background 
metrics reflecting nonzero angular momentum can give rise to axion hair [38] , which in turn 
acts as a source for dilaton hair [39], without the need for U{1) charges. The more general 
case of dilaton and axion hair for Kerr-Newman black holes [40-41] combines all of these 
scenarios . There have also been recent studies of the more physically interesting case of a 
massive dilaton coupled to an Abelian charged black hole [44-45]. It is widely believed (but 
not required) that the dilaton acquires a mass when SUSY is broken: a precisely massless 
dilaton violates the equivalence principle [33,44] , and the dilaton cannot have a mass with 
SUSY intact. Since SUSY is broken at low energies in any event, it seems essential that 
all of the above scenarios be reexamined with a massive dilaton, though the details of the 
SUSY-breaking mechanism and the dilaton potential are not yet well understood. 

The convergence of these separate efforts in black hole physics was inevitable. A natural 
question to ask is whether non-Abelian gauge fields in the low-energy string context lead 
to black holes with primary hair, and if so, whether the sphaleron nature of previous non- 
Abehan solutions is modified enough by the "stringy" scalar fields to yield stable solutions. 
It was the desire to answer these questions, as well as to explore more general black hole 
solutions to what could be the physically relevant theory, which motivated the present work. 
While this paper was being completed, however, we became aware of recent work in Einstein- 
Yang-Mills-Dilaton (EYMD) theory [46-48,49-50] which in part grew out of solutions to the 
Yang-Mills-Dilaton (YMD) system [51-52]. Though these efforts involve strictly massless 
dilatons, some of our numerical results overlap with those of [49] in which the authors 
examine a special case of the more general dilaton coupling 7 explored in [46-47] . We draw 
comparisons to these numerical results wherever appropriate, and discuss the implications 
of this recent body of work for the stabihty of our solutions. 

This paper is organized as follows. In Section 2, we introduce the bosonic part of the 
low-energy heterotic string action, which we take to first order in the inverse string tension 

#2 For reviews of these and related developments, see [42-43] . 
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after compactification to (3 + 1) dimensions. We specify the generic form of dilaton-axion 
potential which arises when SUSY is broken by gaugino condensation, and obtain a simph- 
fied string-inspired theory by requiring spherical symmetry and staticity, and by assuming 
that the characteristic curvature of solutions is small compared to the Planck curvature. 
The spherically symmetric metric and SU (2) connection ansatz are then used to fully spec- 
ify the theory, which is rewritten in terms of dimensionless parameters and variables before 
the general field equations are derived in Section 3. In Section 4, we classify all possible 
non-Abelian solutions to the theory and ignore the embedded Abelian solutions, which cor- 
respond to some of those discussed above but with dilaton and axion masses included. Our 
analysis indicates that only two scenarios can admit solutions: a massive dilaton coupled 
to a single magnetic Yang-Mills degree of freedom, and the full theory of a massive dilaton 
and massive axion coupled to non-Abelian electric and magnetic fields. Though the latter 
theory is numerically intractable, we outline a possible solution scenario before extensively 
analyzing and presenting numerical regular and black hole solutions to the former theory in 
Section 5. In the course of analyzing this theory, which we label EYMD theory, we also note 
the equivalence of scaling arguments for the existence of solutions [53] and a judicious com- 
bination of the field equations. In Section 6, we speculate further on solutions to the most 
general non-Abehan scenario, briefly address the issue of stabihty, and offer our conclusions. 

2. Preliminaries 

Low-Energy String Action 

Our starting point is the bosonic part of the low-energy heterotic string action [54] 



which is expressed in the Einstein frame for a metric with signature (— H — h +). The action 
has been expanded to first order in the inverse string tension a' — jg^^ where g is the 




d X\J —'g 



R 





(2.1) 



gauge couphng for the Yang- Mills (YM) curvature F — dA-\- gA A A and k,^ — SnG. H^yx 
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is the field strength tensor associated with the three-form 

R^dB^ — ^^^L-g^^zv) (2.2) 

where B is the two-form potential in the gravitational supersymmetric multiplet. Q^-^i, and 
are the Lorentz Chern-Simons (LCS) and Yang-Mills Chern-Simons (YMCS) three- 
forms 



^3y — 



uj /\ R uj /\uj /\uj 

3 

1 

A^F--A^A^A 



(2.3) 
(2.4) 



which arise in string theory in order to remove gauge and gravitational anomalies. Here tr 
and Tr denote trace over the suppressed gauge and Lorentz indices, respectively, and the 
normalization for Vlyy is chosen for gauge generators satisfying tr (T*T.^) = —25"^^. R in 
flsL is not the first curvature scalar R^vQ^^ that appears in the action; it is the curvature 
two-form 

R^j,y = du^y u^"" A Uav (2.5) 

where uJanu = {eij.)^Va{eu)b is the spin connection for the tetrad (e^^);,. The other gravita- 
tional scalar appearing in (2.1) is the Gauss-Bonnet (GB) curvature combination 

R^ = Rf^ypaR^^'P'' - m^uR^^ + (2.6) 

which also helps to cancel anomalies and is second-order in derivatives of g^jj. The dilaton 
field D couples to other fields through exponentials with coupling strength 7 and has a self- 
interaction V whose form will be specified below. The normalization of V has been chosen 
to accomodate a choice of coupling and a field rescaling: we take 7 = \/2k and define the 
dimensionless dilaton field = kD/\/2, so that the dilaton kinetic and potential terms 
assume the form 

1 

2^ 

It is important to note that 7 is the only coupling parameter that we fix in our analysis; all 
others (including the k?' factor now appearing in front of V) will be absorbed in the definition 
of other dimensionless fields and parameters. 



-2dn<pd^<p - Ak^V (0, s) (2.7) 
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The field s appearing in V is the dimensionless pseudoscalar Kalb-Ramond axion, the 
only truly dynamical mode of the three-form field which we introduce via 



(2.8) 



With this relation and the dual of the Bianchi identity 



dH= — [Tr {RAR)- ghi {F A F)] 



(2.9) 



which follows from (2.2), we can express the three-form field strength as a sum of axion 
kinetic and topological current contributions 



2k' 



(2.10) 



where ^ {*Hf = (a'/S) [V^ C^sl)^ - g^V ^ {*^wT] is comprised of the four-divergence 
of two topological currents 



V^(*Jl3y)'' = -2Vptr 



(2.11) 
(2.12) 



which can also be expressed — Tr {^R^iiR'^^^ /2 and —tr [f^vF'^^^ /2, respectively. With 
this replacement, the action becomes 



2«2 



R - ^e^^d^sd^'s - 2d^<pd^<i) - Ak^V (0, s) 



+ 



a' (\ 



4 V4 



-s e 



pXjiv 



^aPpX^lJLV + ^9 ^pX^fJ-i^ 



e 

2 



-2(f> 



(2.13) 



The above equation is a general expression for the low-energy heterotic string action to first 
order in the inverse string tension for 7 = y/2K. We now briefly examine some features of 
this form of the action while arriving at some useful simpliflcations. 
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Note that when V = 0, all sources for the dilaton and axion fields are 0{a'), so the fields 
themselves are first order in a'. Furthermore, the sources are comprised of gauge field and 
higher- derivative curvature combinations on an equal footing. Even for fixed s, for which the 
topological current terms arising from the YMCS and LCS three-forms contribute nothing 
to the equations of motion, we must include the GB term if we are to account for the 
gauge field strength. Prom this perspective, the Reissner-Nordstrom solution (corresponding 
here to the Abelian sector of some non-Abelian field strength and fixed dilaton field) should 
for example be viewed as an 0{a') correction to the Schwarzschild solution, subject to GB 
curvature corrections at the same order [41]. The gravitational effect of such curvature terms 
for fixed (f) has been examined in [55-57], although in d = 4 the GB contribution enters purely 
as a boundary term and can be ignored. The inclusion of the dilaton, however, introduces 
an source term in the dilaton field equation even in d = 4, and such scenarios have also 
been studied [58-62] . Several authors have neglected the GB curvature contribution in their 
investigations of the dilaton while consistently keeping the gauge field strength source [31-33] 
by considering solutions whose mass scale is large compared to the Planck mass. In some 
circumstances, such as the extremal limit of charged dilaton black holes [31,63-64], one can 
satisfy the mass scale assumption but introduce a different inconsistency: in this regime, a' 
is necessarily large, so the dropping of higher order terms in the effective string action is 
no longer justified [32]. Mindful of these concerns, we neglect the R^ term in the action by 
assuming that the mass of solutions is large relative to the Planck scale, but that a' is small 
enough for (2.13) to remain reliable. 

Por a dynamical axion s with or without the dilaton, we again encounter higher-order 
curvature and gauge field source terms. Por spacetimes with rotation, the LCS combination 
gives nontrivial contributions to the dilaton and axion equations of motion, and analytical 
solutions for dilaton and axion hair outside Kerr [38-39] and Kerr-Newman [40-41] black 

holes have been obtained. Por the 4-d Schwarzschild spacetime, spacetimes related by a 
conformal transformation, or any 4-d spacetime with a maximally symmetric 2-d subspace, 
the LCS three-form either vanishes or is exact [65]. Thus for the static, spherically symmetric 
spacetime we investigate below, the remaining O (^a'R^) term in the action can be ignored. 

The inclusion of the potential V introduces an additional mass scale into the problem, 
so that the dilaton and axion need not be 0{a'). None of the other preceeding observations 
are qualitatively altered by its inclusion, but in choosing a potential our discussion must 
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move from generic features of heterotic string theory (in static, spherically symmetric (3 + 
1) dimensional spacetime) to a more specific model. Before doing so, we summarize the 
simplifications outlined above by rewriting the action 



R - ^e'^'^d^sd'^s - 2d^(j)d^'(j) - Ak^V (0, , 



a 
T 



(2.14) 



where we have used the fact that F^^F^^ is a 4-divergence ( eq.(2.12) ) in a topologically 
trivial spacetime to recast the axion-gauge field coupling in a more convenient form. 

Dilaton-Axion Potential 

We choose a potential of the form which arises when supersymmetry is broken by gaugino 
condensation in the hidden sector of the theory [66] 

^2 



y(0,.) = /^ 



a 



1 



(a— aoo) 



+ 11 



2,i^±^e-(«-"-)/2cos( A, 
(«oo + l) V2&0 



(2.15) 



where a = 3exp(— 20)/6o, ccqo corresponds to the dilaton field at a potential minimum, and 
bo is determined by the one- loop /3-function of Q, the subgroup of the hidden sector gauge 
group which precipitates supersymmetry breaking. We take Q to be the entire hidden sector 
gauge group that arises in these scenarios, £"3, for which bo = 90/(167r2). The parameter n is 
a scale related to the vacuum expectation values of the gaugino pair xx the three-form 
Hmnp-i where m, n, and p are indices on the internal compact manifold K only; we treat it 
here as a free parameter. With the axion field set to its vacuum value, Sqo = &o(4vrn)/3 for 
integer n, this potential has been used by some authors to investigate infiation and cosmology 
in the context of superstring theories (see e.g. [62, 67-69]). A plot of V lor s — Sqo and 
0OO = is shown in fig. la; it has a minimum at = 0oo and achieves a local maximum at 
(f) < before y — > for — > —00. 

Metric 

We parametrize the metric for a static, spherically symmetric spacetime as 

ds^ = (r) df + (r) dr'^ + {dO'^ + sin^ Qdip^) , (2.16) 

where R (r) = (1 — 2Gm (r) /r)~^^^, m (r) is the total mass-energy within the radius r, and 
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we have set c = 1. To describe black hole solutions, we define 6 = — In (i?/T) and rewrite 
(2.16) as 

1 _ dr' + {dO' + sin^ ed^') . (2.17) 

Regularity at the origin requires T (0) < cxo and R' (0) , T' (0) = , while regularity at the 
event horizon at r = r/j is satisfied by m (r/J = r/i/2G' and 6 (r/j) < oo. Wc also impose the 
condition of asymptotic fiatness, which implies i? (r) , T (r) — 1 as r — >■ oo or, equivalently, 
R{r) ^ 1 and S (r) 0. By exploiting the freedom to rescale the time coordinate, however, 
we can make the boundary conditions more suitable for integrating the Einstein equations. 
Rather than requiring Tq = T (0) and 5q = S (r/j) as initial conditions, consider rescaling t 
such that 

it:(o) = i, T(o) = i 

(2.18) 

The condition of asymptotic flatness then translates into ^(oo) = 1/Tq and S {oo) — —Sq, 
and we can determine Tq or 5o from the behavior of a solution as r — > oo. 

Gauge Connection and YM Curvature 

In this paper, we investigate the simplest choice of non-Abelian gauge group, SU{2). 
The most general spherically symmetric SU (2) connection can be written in the form [70] 

A = - [oTr dt + hTr dr + idre - (1 + c) f^) d6l + ((1 + c) f^i + dr^) sin Od^p] , (2.19) 
9 

where g is the gauge coupling and (t>, tq-,t^p) is the antihermitian su{2) basis projected along 
the polar coordinate directions: ?> = r • r, etc., and the matrices satisfy [Ta,r6] = eah(^c (we 
deviate from the gauge generator normalization used above for this section only). From 
(2.19) we note that c and d are in general dimensionless functions of r and t, while a and h 
have dimension The connection has a residual gauge freedom under transformations 

of the form U — exp [(5 (r, t) tv], where (5 (r, t) is an arbitrary real function, which we can use 
to set 6 = while preserving the form of the connection. Because the held equation for b 
becomes a potentially useful constraint equation after gauge fixing, we leave h nonzero for 
now and demand that the component functions of A depend on r only. 
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Following 
scalar form, 



, we express the c and d degrees of freedom in the connection in complex 



c (r) — id (r) = / (r) exp i(5 (r) , 



(2.20) 



which will make the non-Abelian character of the system more transparent. The YM curva- 
ture F — dA + gA A A for the connection (2.19) is then 
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af-dt +{b- (3') -^dr 



1 

+- 
9 



rR 
f 

— dr 
rR 



+- 



(cos Pro — sin /3t^) d9 + (sin Ptq + cos /3f^) d(p 
(sin + cos Pt^p) dO + (cos Ptq — sin /Sr^) dip 



-a'—Trdt A dr + ^ (f - l) Trd6 A dip 
R ^/ V / 



which we have expressed in a convenient orthonormal tetrad basis, 



dt = —dt, dr = Rdr, d6 = rd9, dip = rsvaOdip. 



(2.21) 



(2.22) 



Note that the dependence of F and A on the gauge coupling indicates that g only enters the 
simplified form of the action (2.14) through the inverse string tension a' . 

3. General Field Equations 

With the choice of ansatz (2.19) for the gauge connection and (2.15) for the dilaton- 
axion potential, we are in a position to express the action in terms of the axion, dilaton, 
metric, and gauge degrees of freedom. Before proceeding, we relax the condition c = 1 and 
examine the dimcnsionful quantities in the action in order to cast our equations of motion 
in dimensionless form. 

Noting that [g] = [r][M]-V2[L]-3/2, [^,2] ^ [r]-i[M] V2[i,]-i/2^ and that k = V^/c^ 
has dimensions [n] — [T][M]~-^/^[L]~-^/^, we observe that the parameters appearing exphcitly 
in the action have dimensions 



[a'] = [n'/g^] = [L]\ [«:V] = [L]-\ 



(3.1) 



and that from (2.13) all terms in 2k C have dimensions [L] ^. Thus, if we factor {a' /-i) 
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out of 2k^C, define 'g = 2/\/a', and define the dimensionless quantities 



= i^ii^jg, r = gr, a = a/g, h = b/g, 



(3.2) 



then ^/=g(2«;2£) 

can be written purely in terms of dimensionless fields and parameters. To 
do this explicitly, we also define a dimensionless mass-energy function m based on the metric 
(2.17): 



i?2 (r) 



1 - 



1 - 



2m 



for fh (f) = gG m {gr) jc? . 



(3.3) 



Expressing the curvature scalar R^j^^g^^ in terms of the metric functions i?(m, r) and T'(r), we 
find the following expression for the gravitational and matter action of our static, spherically 
symmetric system: 



Sg =^ I dtdr 



2\B? dr\T 



gG 



dtdr -m4 ( e-^) (3.4) 
^2 ' 2 i?2 +r'V{<f>,s)y^-s'a{f-l) (3.5) 



Sm / dtdr 

gG J 

+ e-2<^ ( 



1 {ra'f 

2 i?2 



1 



where the prime denotes derivative with respect to r and V = k^V. The solutions to the 
dimensionless field equations obtained from this action give us solutions for any ^ > ( or 
< a' < oo) through scahng relations. 



«? {r) = ga (gr) , (r) 



gG 



(gr) , (r) = T {gr) , 



(3.6) 



where denotes any of the functions {s, 0, /, i?, T, (5} and we have ignored 6 since it will 
be eliminated by gauge fixing. Hence the radial structure of solutions for a given value of 'g is 
the same as that obtained from (3.5), but it occurs at a physical radius r = r/^with physical 
scales given by (3.6) and //^ — 'gji'^ / n. For notational simplicity throughout the remainder 
of the paper, we drop the carets on dimensionless quantities with the understanding that 
everything is now dimensionless unless otherwise specified. 
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By varying (3.5) with respect to the fields, we obtain the dimensionless, static field 
equations 



dr\RT J d(f) ^ T 2^ RT 



(3.7) 



+ 2e 



-2(j) 



i?2 



+ f iP' - b) 



2 i?2 + 



(^e-2^|r2a'^ - 2e-^'^faRT - s' {l - f) = 



R 
T 







dr 
d_ 
dr 



dr \ RT 



-/ 



R 



i?2 



f + T'a'f\--s'af^O 



dr \ RT 



f [P'-b] ^o, 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



and the constraint equation 



'rt 



./2 - 6) = 0. 



(3.12) 



Even without using the remaining gauge freedom to set 6 = 0, we find that /2 (/?' — b) 
disappears from the field equations. The constraint equation with gauge fixing does give us 
additional information, however: it implies that our gauge choice eliminates an additional 

degree of freedom (c and d are related by a multiplicative constant) when we insist / ^ 0, a 
criterion for non-Abclian solutions. In the Abclian sector / = 0, it and the f3 field equation 
indicate that f3' is an arbitrary function of radius, which refiects the fact that c = d = in 
(2.20) entirely eliminates the the need for the complex scalar phase. In either case, we can 
now eliminate p {(5' — b) from our analysis. 

To obtain the Einstein equations, we can either utilize the energy-momentum tensor. 



-2(f) 



+ 



+ 2d^(f)du(f) - 



-e^'^dpsdPs + dp(t)dP(t> + 2V (0, s) 



(3.13) 



or use the explicit dependence of (3.4) and (3.5) on some pair of independent metric functions 
(m and 5 will do as well as R and T) to derive the gravitational equations directly. By 
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whatever route, we find the {tt) and (rr) Einstein equations can be expressed in the form 

,2 



m' = 6-2-^ 



1 



2m\ 



2m ^ 



+ 11-^ 

r 



(3.14) 



2m\ T' 

1 = e 

r 7 T 



2.;!- 



- 1- 



r / 

2m 



2r2 



e40(rs')Vi(r</.')' 



2 ™ 



+ V (0, s) - 



(3.15) 



For black hole solutions, we replace the T' equation with 

-2 



r 



-2(j) 



/2+(i_^y ^25^2^2 



(3.16) 



which indicates that 6 monotonically decreases with radius. From this fact and the boundary 
condition 6 (oo) = 0, we can conclude that 6 > and thus T {r) > R{r) > 1. Another metric 
relation useful for predicting solution properties is 

,2- 



d 

dr 









2e-2<^ 


'r 









, (1-/2)' 



1-— r+ 



2r2 



+ 2e" 



-2^ 



- 2V (0, s) 
RT 



R 
T 



(3.17) 



which may be obtained by combining the field equations. In preparing the equations for 
integration, we will also require the relation 

>2- 



RdrXRT 



2m — e 



2V (<l>, s) 



(3.18) 

Finally, in classifying the possible solutions of this system of equations, it is useful to rewrite 
(3.9) as 



1 d 
2dr 



e-^-*l iraf + 2e-^fVRT + s' (l - f) a. (3.19) 
R 



The positivity properties of the right-hand side of this equation, coupled with boundary 
conditions, can be applied to establish no-dyon results in the non-Abelian sector analogous 
to those for EYM [11] and EYMH [23] theories. 
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4. Taxonomy of Non- Abelian Solutions 



In previous studies of SU{2) gauge theories coupled to Einstein gravity, many authors 
have noted that / = in the gauge connection ansatz gives a theory with Coulombic (i.e. 
U{1)) magnetic and electric charges. As discussed in the introduction, a great deal of work 
has recently been done on the corresponding Abelian sector of the low-energy heterotic string 
action. Though the addition of the SUSY-breaking potential (2.15) could provide interesting 
new features for some Abelian sector solutions (see, for example, [44-45]), we choose not to 
explore them here. Our focus is the non-Abehan sector of the theory (2.14), and in this 
section we consider the various possibilities for static, spherically symmetric solutions with 

The most general class of possible non- Abelian solutions to the field equations corre- 
sponds to {/, a, s, 0} behaving as nontrivial functions of radius. By studying the asymp- 
totic behavior of the field equations, we find that the condition / ^ requires (oo) = 1, 
which further imphes a (oo) = via (3.9)-(3.10). The decay of / and a toward these asymp- 
totic values is also too rapid to give a net electric or magnetic charge, so such solutions are 
not dyons in the usual sense but do possess both electric and magnetic fields. The asymp- 
totic values of and s are fixed by V; as we observed in the discussion of the dilaton-axion 
potential, we are free to choose 0oo but the axion must assume the value Sqo = 6o(47rn)/3 
for some integer n. At the origin or horizon of a black hole solution, we must choose 
{/" (0) , a' (0) , s (0) , (f) (0)} and {/ (r/^) , a' {rh) , s {vh) , (t> (r/i)}, respectively, to integrate 
the equations of motion; in either case finite energy density restrictions force the initial 
value of a to vanish. Exploring a system with four integration parameters is impractical 
numerically, but we can demonstrate analytically that solutions are not forbidden to exist, 
and we can anticipate some of the properties of potential solutions. The analysis of the 
general case is best done, however, after the simpler possibilities are surveyed. To describe 
the less general non-Abelian solutions, we reduce the number of integration parameters by 
ansatzing each of the fields {/, a, s, 0} to an appropriate constant value in turn. We will 
find that only the a = constant and s = constant cases admit nontrivial solutions, though 
the = constant ansatz provides much insight into the characteristics of the most general 
non-Abelian solutions. 

f = constant 
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The asymptotic behavior of the field equations restrict the acceptable values to / = 1, 
and the alternate version of the a equation (3.19) can then be used to establish a = 0. The 
proof relies on the fact that the right-hand side of (3.19) is nonnegative for = 1, so that 

is strictly increasing beyond the initial radius for nontrivial solutions. The boundary 
conditions a (0) — a (r/j) = a (oo) = 0, however, require that decrease toward = 
asymptotically. It then follows that the only solution consistent with the a field equation 
and the boundary conditions is the trivial solution a = 0, and we are left with a theory 
of two gravitating scalar fields. Through the metric relation (3.17) and the same line of 
reasoning used to establish a = 0, we can now show that only trivial solutions follow from 
/ = constant. Since the right-hand side of that relation is never positive, nontrivial solutions 
obey the condition T' > beyond the initial radius. Because we also noted in the discussion 
of (3.16) that T (r) > T (oo) = 1, the only way to reconcile monotonically increasing T with 
the boundary conditions is to require (0, s) = 0: only the trivial solution = 0oo, s = Sqo 
is compatible with both the boundary conditions and the field equations. 

a = constant or s = constant: EYMD Theory 

By setting either a or s to a constant value, we arrive at a theory with one gauge degree of 
freedom coupled to a massive dilaton and Einstein gravity, which we denote Einstein- Yang- 
Mills-Dilaton (EYMD) theory. In the s = constant case, the potential requires s = Sqo, 
and the a equation (3.19) with s' — subsequently yields a = when we employ the 
reasoning introduced in the / = constant discussion. Alternately, if we take a = constant 
the boundary conditions require a = 0, and the field equations then imply s = Sqo for ^ 1. 
Note that the right-hand side of the metric relation (3.17) again involves a —V contribution, 
but now there are positive /-dependent terms which make possible the decrease of T (r) 
toward T (oo) = 1 required for nontrivial solutions. To obtain regular or black hole solutions 
to EYMD theory, we must fix ipoo and choose the integration parameters {/" (0) , (0)} or 
{/ (f^h) ) (^/i)}i respectively, such that the fields match the appropriate boundary conditions 
upon integration. We provide a detailed discussion of this procedure, which is often referred 
to as a two-parameter "shooting" procedure, and obtain numerical solutions in the next 
section. 

It is interesting to note that in the absence of the dilaton-axion potential, solutions to 
EYMD theory are the most general non-Abelian solutions to the full field equations: dyonic 
non-Abelian solutions and nontrivial axion solutions are prohibited when V = 0. The s field 
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equation (3.8) in this case implies 

RT 



' = 4^e-^^ [«(!-/')+ ^^']. (4-1) 



s 

r 



where Cg' is a constant which acts asymptotically like a Coulombic charge, giving a contri- 
bution to the total mass-energy proportional to c^z/r^ at large radius. Though nonzero Cg' 
appears acceptable at large r, it violates regularity of the metric at the origin for regular 
solutions and at the event horizon for black hole solutions, since m' (0) and 5' (rh) diverge 
as {rs') diverges. When (4.1) with Cs' = is substituted into the alternate form of the a 
equation (3.19), the right-hand side is nonnegative and a = by the arguments outlined 
above. It then follows from (4.1) that s = constant, the value of which is no longer fixed by 
the minimum of V, and the general theory reduces to EYMD theory for a massless dilaton. 
A closer look at the field equations, however, reveals that a = and s = constant also follows 
from the weaker condition (dV/ds) = 0: the no-dyon result is a consequence of having an 
axion-independent potential, rather than no potential at all. The possibility of non-Abelian 
solutions other than the EYMD class relies crucially on the presence of a massive (or at least 
self-interacting) axion field in the theory. 

(f) = constant 

At first glance, setting equal to its asymptotic value 0oo appears to reduce the number 
of independent matter fields from four to three. The exponential coupling of the dilaton 
to the gauge field and axion kinetic terms in (3.7) instead make the dilaton field equation 
a nontrivial constraint which must be satisfied by the remaining fields {/, a, s}. Using 
this constraint and some field redefinitions, it is possible to cast the resulting theory in a 
form which requires only two independent integration parameters. Though the constraint 
equation appears to make the theory numerically tractable, it is also responsible for the 
nonexistence of solutions: by differentiating with respect to r, we find that the constraint is 
incompatible with the = constant field equations. Because this theory shares important 
features with the most general non-Abelian case, we develop the two-parameter formulation 
and show explicitly that solutions cannot exist. 

The prospect of black hole solutions with nontrivial gauge degree of freedom a motivates 
us to define A = aRT: based on the expression (3.16) for 5', a must vanish at least as fast 
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as 1/ RT near the horizon for the a? term not to diverge and violate regularity requirements. 
It is also convenient to introduce £ = e^a' in lieu of A! ., to which £ is related by 



(4.2) 



In terms of these new variables, the self-interaction and kinetic contributions of a to the action 
interchange roles; the factor which appears with (a')^ is absorbed by the definition of 
and reappears in the f'^o? term: 



2^ 
i?2 



(4.3) 



This combination also appears in the Einstein equations (3.14)-(3.15) and the metric re- 
lation (3.17), which are unchangedapart from the substitution of (4.3) and = 0oo- The 
interchange of kinetic and potential roles is sensible when we consider the expression for 5', 



-20O 



(4.4) 



which now involves only kinetic terms and is no longer explicitly dependent on the metric 
functions. Because the relation between £ and A' depends explicitly on 5', it is convenient 
to use (4.4) in place of the Einstein equation (3.15) when considering regular as well as black 
hole solutions. 

Through the metric derivative equation (3.18), which becomes 



2 2r2 



^(0oo,s)r2 



(4.5) 



for this case, we can also express the axion and gauge field equations in a form independent 
of either T or 5, 



2r 2 5± 



s' - 4e-^'^^i 



4e" 



dV{<P,s) 
ds 



^ (1 - f) - 2ff^ 



(4.6) 



= 
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-"^S' + 2r8 - 2fA - e^^^s (l - /^) = 



i?2 



2 (5 

r e 



d f e 



dr V /22 



(4.7) 
(4.8) 



The potential for (j) = (p^Q may be written in the form 

,1/3 

V (0OO, s) = 2jj,'^ ( 1 - cos ' 



4/sin2 



.2 V^'o 

which is proportional to the derivative with respect to at (poo 



1 /3 



4 \bo 



(4.9) 



apo {aoo + 1) + 2 
(aoo + 1) 



V{^oo,s), (4.10) 



where ctoo = 3 exp(— 20oo)/&o and we have introduced A^^ to denote the constant of propor- 
tionality. This property of the potential allows us to write the dilaton constraint equation 
in the virial-like form 



2r' 



2^^ ^ 



0. 



(4.11) 

We can use this constraint to help integrate the system comprised of the m' and 5' equations, 
eqns.(4.6)-(4.8), and the E-Al relation (4.2): it reduces the number of integration parameters 
to two and provides a check as we integrate the system. 

There are reasons to expect nontrivial solutions to this theory. We can formally integrate 
(4.6) to obtain an expression analogous to (4.1), 

\2 



Moo (^^0^ 



eh' 



dr r 



ds 



(4.12) 



which upon substitution into the alternate form of the a equation (3.19) gives us 



d_ 
dr 



2e 



= 2e" 
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2^^ ^ 7?2 

2 ^2- 



, 1 , 



A 



1^2^ ^ + i?2 j + 8 



Mo 



\rs 



(4.13) 



— s 



OS 



In the two previous cases considered, we used the manifest positivity of the right-hand side of 
this equation to help establish A — £ = Q and s = constant as the only acceptable solution. 



19 



Now we must contend with the final term in (4.13), the sign and magnitude of which depend 
on the details of the potential. For the pure axion potential depicted in fig. lb, we can 
imagine a solution scenario in which s increases from some initial value 47™(6o/3) < sq < 
47r(n+ l/2)(6o/3), moves over the potential maximum at Smax — 47r(n+ l/2)(6o/3), and 
comes to rest at the closest degenerate minimum Soo = 47r(n+l)(6o/3). With s' > along the 
entire solution trajectory, the sign of the final term in (4.13) is completely determined by the 
integral factor. For Smax < s < Sqq, the integrand is negative-definite. For sq < s < Smax, the 
sign of dV/ds changes but the integrated contribution from r > r{smax) initially dominates. 
Since the monotonically increasing factor r^e~^ in the integrand weights the r > r{smax) 
contribution more heavily, and \dV/ds\ is symmetric about s = Smax, we can in fact conclude 



OS 



< 



(4.14) 



for the entire trajectory, and the right-hand side of (4.13) is not positive definite. In order that 
a be nontrivial, however, the magnitude of the contribution (4.14) must be large enough that 
the integral of the right-hand side of (4.13) be negative as r — > oo, which from the original 
equation (3.19) is required for the asymptotic decrease of toward zero. From (4.12), we 
must therefore have 



4 i?2 



{rs'Y 



e's' 



f 

Joe 



dr r 



ds 



(4.15) 



at least as r ^ oo, from which we conclude ^ (l — Z^) < asymptotically and A < if f'^ 
approaches (oo) = 1 from below. We can infer some general features of by rewriting 
the / field equation as 



2 dr 



i?2 



if)' 
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,f_ 
i?2 



— e 



(4-16) 



From this expression, we find that solutions must obey the inequality 



(4.17) 



when ^ < or ^ > e'^'^°°s', since outside this region the right-hand side of (4.16) is positive- 
definite and /2 cannot approach unity asymptotically. For < ^ < e'^'^°°s', on the other 
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hand, can exceed the bound in (4.17) but must then approach (oo) monotonically, 
so the incquahty simphfics to < 1. Further imphcations follow from the original form 
of the / equation (4.8), which indicates that //" < when /' = and (4.17) is satisfied 
regardless of the value of A: solutions can exhibit oscillations about / = at finite radius, but 
they cannot possess turning points once the inequality bound is exceeded. Combining these 
observations, we can construct a simple picture for the s' > trajectory that is consistent 
with the field equations: ^ < 0, with a single extremum and increase toward A (oo) = as 
the negative contribution (4.14) begins to dominate in (4.13), and ^ 1, with oscillations 
and nodes in / possible. Had we begun with 47r(n + l/2)(6o/3) < so < 47r(n + l)(feo/3) and 
demanded that s decrease monotonically toward Soo = 47™(6o/3), we would have found an 
identical picture with A —A. Thus, in the presence of a sufficiently massive axion field, 
non-Abelian dyonic solutions for the 4> = constant theory appear possible. 

Unfortunately, the constraint (4.11) that follows from the = constant ansatz is not 
consistent with the rest of the field equations. It is not obvious from inspection whether 
solutions to (4.6)-(4.8), (4.2), and the Einstein equations satisfy (4.11) for some choice of 

the parameters {/j, (f)oo}- Because this set of equations without the constraint is sufficient to 
determine solutions, we should be able to verify (4.11) by utilizing the entire set. When we 
take the derivative of (4.11) with respect to r and use the other field equations to eliminate 
higher derivatives, we obtain the relation 



2e 



\ 



(f)'^ + 2Aef] - 3/ ((1 - /^)' + (1 _ 

J 



f 

' (4.18) 

which should be satisfied along with (4.11) for all r. It appears that no choice of {/x, 0oo} 
or further simplification with the field equations can make (4.18) an identity for nontrivial 
gauge and axion fields, and we conclude that the (f) = constant case admits only trivial 
solutions. 

General Non-Abelian Solutions 
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The failure of the = constant system to admit nontrivial solutions is closely tied to 
the constraint created by a nondynamic dilaton field. How does the situation change when 



we relax the restriction that 6 = 



Exciting the dilaton degree of freedom restores the original form of the potential (2.15), 
but the expansion of the potential in powers of = — can be written in the form 



V 



(a^ + ttoo + 2) (a^ + + 6aoo + 4) ^ 
1 + ^ — ; —7^ (p+- o -(p +U 



(aoo + 1 



2 (aoo + ly 
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(4.19) 

Thus dV/ds deviates from the purely axionic form discussed above, but only by a 0- 
dependent scale factor. This scale factor breaks the symmetry of \dV/ds\ about s — Smax, but 
it is of order unity when is near the minimum at 0oo) so the reasoning behind (4.14) for the 
s' > solution scenario remains intact when the expansion (4.19) is valid. Relaxing = 0oo 

also gives dynamical exponential couplings in the field equations, but neither these changes 
nor the scale factor seem to qualitatively alter our solution discussion in this regime. The 
dilaton kinetic term is restored in 5' , 



-2, [;.2 + ;2^2] + 1^4,(^^.)2^1^^^,y 



(4.20) 



and similarly in the expressions (3.14)-(3.15) for m' and T', while the axion and gauge field 
equations acquire additional terms linear in 0': 



s" + 
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2r + 4rV) ,d 

i?2 ^ dr\B? 



dV{<t>.s) 
ds 



(4.21) 
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-'^8' + (2r - 2r20') S - 2fA - e^^s' (l - Z^) = 
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2r 



2a> 



f' + 



(4.22) 
/ = 0. (4.23) 



The equations (4.12)-(4.13) and (4.16), which were derived from the axion and gauge field 
equations, retain their form, and the inequalities (4.14)-(4.15) and (4.17) which followed 
from them are still valid if we take 0oo — 4> {r). 
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Though the solution scenario explored in connection with (f) = 0oo appears promising 
when the scale factor in (4.19) is convergent, the general non-Abelian case still requires an 
analysis of the dilaton equation 



which we have rewritten in a form analogous to (4.16). Since the inequalities derived above 
provide no information about the relative magnitudes of the terms on the right-hand side 
of (4.24), it is difficult to extract an inequality restricting the behavior of 0, but we can 
make some observations. Since dV/d(t) < for < and near the minimum, the potential 
contribution and possibly the gauge field contribution to the right-hand side of (4.24) are 
positive, which tends to drive away from the vacuum. For > near the minimum, these 
signs reverse and make a monotonic decrease of toward = more likely. Though our 
preceeding analysis does not require such monotonic behavior for the dilaton, it is certainly 
consistent with the solution scenario and provides a potentially viable alternative to the 
(f) = ansatz. 

The complexity of the full string-inspired theory is well reflected in our inability to 
proceed any further analytically. It appears that the question of existence of non-Abelian 
solutions with nontrivial dilaton and axion fields can only be settled by further numerical 
study. Though we have not actually obtained numerical solutions in the general case, we have 
considered a strategy for simplifying this four-parameter problem. In past (two-parameter) 
shooting problems, we have found the method of "shooting to a fitting point" a convincing 
way to confirm our numerical results. In this method, shooting parameters at both fixed 
points (r = 0; r/j and r = oo) are chosen so that the field equations can be integrated toward 
a common midpoint, at which the two trajectories join smoothly if the parameters corre- 
spond to a solution. Though this procedure requires more than double the original number 
of shooting parameters, the deviation of the trajectories at the midpoint can be used to 
choose new shooting parameters via a multi-dimensional Newton-Raphson algorithm. If the 
initial shooting parameters are reasonably close to a solution, this procedure is moderately 
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(4.24) 
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successful at converging on the solution, but it is difficult to match the trajectories with an 
error comparable to the global tolerance of the integration routine. Once the neighborhood 
of a solution is determined, this appears to be a more promising approach to finding the 
solution than the method employed in the next section, at least for three or more shooting 
parameters. We hope to utilize this procedure in the future to find solutions to the full 
string-inspired theory. 

By closely examining some simplifying ansatzes, we have been able to narrow the possible 
solution classes to two: a massive dilaton coupled to a single magnetic Yang-Mills degree 
of freedom, which we denote EYMD theory, and the gauge field coupled to both massive 
dilaton and massive axion fields. Though we could only speculate about solution scenarios 
for the latter theory, in the next section we present and analyze numerical solutions to the 
former. 



5. EYMD Theory: Regular and Black Hole Solutions 

EYMD Equations 

As we observed above, the ansatzes = or a' = in the non-Abelian sector give a = 
and s = Soo for both regular and black hole solutions. To obtain numerical solutions to the 
resulting EYMD system, it is instructive to reexpress the dilaton as (f){r) = h{r) /r + 0oo, 
where h/r, which we also denote 0, is the deviation of from its vacuum value. With this 
change and the simplifications a = and s — Sqo — ^'o(47m)/3, the remaining gauge field 
equation becomes 



dr \RT 



and the dilaton equation assumes the form 



T 



RT 



(5.1) 



where V = dV/d(t)\s=s^. The Einstein equations (3.14)-(3.15) for EYMD theory are 



r / 



(5.3) 
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with the auxihary equation 
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r 



(5.5) 



replacing the T' equation for black hole solutions. To simplify the integration of the equations 
of motion, we use the metric derivative relation (3.18) to express the gauge field and dilaton 
equations in a form independent of T or 5: 
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(5.7) 



Though it is the /i-dependent form of the field equations we choose to integrate, we pass 
freely between 0, (f), and h in our analysis when describing the dilaton. 

Analytical Features and Boundary Conditions 

We can anticipate the general features of the solutions and boundary conditions from 
the field equations.The gauge field equation may be rewritten 



1 d 
2dr 
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RT 
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(5.8) 



Since the right-hand side of this equation is manifestly positive for > 1, the only nontrivial 
solutions having finite must satisfy < 1. By expanding the left-hand side of the 
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equation, we can see that solutions satisfy //" < for /' = and < 1, which is 
characteristic of oscillations about / = 0. As in past studies with non-Abelian gauge fields 
coupled to gravity (see e.g. [7-10]), we expect that solutions will be classifiable by the number 
of nodes which occur as / oscillates. 

Prom the / field equation, we note that the boundary conditions for the gauge field 
include |/| = 1 and / = as r — > oo. The asymptotic behavior of the field equations reveal 
that the latter condition implies / = 0, which for black hole solutions corresponds to an 
Abelian gauge field carrying purely magnetic charge. Without the dilaton, this case just 
reduces to the Reissner-Nordstrom solution, while with the dilaton we recover the type of 
theory studied in [44-45]. The / = case for regular solutions is forbidden by boundary 
conditions at the origin: regularity of the metric only allows the possibility |/ (0)| = 1. If we 
attempt to set |/| = 1, corresponding to the theory of a massive dilaton coupled to Einstein 
gravity, we find that both regular and black hole solutions are forbidden: according to (3.17), 
T must be monotonically increasing, but boundary conditions provide the incompatible 
restrictions T (0) > T" (oo) and T (rh) > T{oo). We conclude that for fundamentally non- 
Abehan solutions to EYMD theory, we must have \f (oo)| = | / (0)| = 1 for regular solutions, 
1/ (oo)| = 1 and / (r/j) unspecified for black hole solutions, and (r) < 1 for all solutions. 

To better understand the expected behavior of 0, we rewrite the dilaton equation 
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where = (h/r) is the deviation of from the asymptotic value fixed by the minimum of V. 
For < near the the minimum, l^' < and the right-hand side of (5.9) is positive-definite, 
which implies that is strictly increasing. Although V' changes sign before it vanishes as 
— > — oo (cf. fig. 1), the gauge field contribution to (5.9) is exponentially amplified in the 
same limit, so the details of the potential should not alter the conclusion that the region 
< is forbidden. We can also predict that solutions exhibit a monotonic decrease from 
> to the vacuum = 0. This is unambiguous in the /j,'^ — case, for which (5.2) indicates 
0' is strictly negative, but nonzero V' introduces the possibility of turning points. Carrying 
out the derivatives in the left-hand side of (5.9) demonstrates that for 0' = 0, 0" < if the 
gauge field contribution dominates V', while 0" > iiV' is dominant. It follows that (j)' — 
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can only occur at a plateau in the former case, but for V' dominant a turning point in is 
possible. Since the gauge field contribution is exponentially suppressed and V' grows even 
larger as (p increases, it appears that a turning point is a precursor to diverging and cannot 
be a solution feature. To summarize, we expect all solutions for which (oo) is finite to be 
characterized by (f) rolling monotonically to zero from above, and |/| approaching unity after 
crossing at least once through f — 0. 

In addition to providing a sketch of the behavior of the fields for regular and black hole 
solutions, the preceeding analysis provides insight into the behavior of the fields when we 
are close to a solution. To qualify what we mean by "close", we introduce some details of 
the approach to solving this system numerically. For regular solutions, finite energy density 
{k^ /'g)TQQ and regularity of the metric at the origin give the following behavior as r — > 0: 

/(r) = -l + ^r(0)r2 + O(r4) (5.10) 

h{r)^h' (<d)r + 0{r^) (5.11) 

2m{r) = 0{r^) (5.12) 

Inr(r) = C>(r2), (5.13) 

where /" (0) and h' (0) = (0) are parameters we must adjust to "shoot" a solution to 
match the asymptotic conditions |/ (oo)| = 1 and (oo) = 0, and we have taken / (0) = — 1 
and used the rescaled initial condition T (0) = 1 introduced in (2.18). All the terms not 
shown explicitly in (5.10)-(5.13) depend only on the two shooting parameters and the free 
parameters /x^ and 0oo associated with V . Similarly, for black hole solutions we can use the 
metric condition 1/i?^ (r/j) = for m (r/j) = r/i/2 to expand near the horizon: 



f{r) = f {rn) + f {ru) (r - vu) + O {r - ryf , (5.14) 

h{T) = h (rft) + h' {n) (r - r/0 + O {r - rnf , (5.15) 

m (r) = rh/2 + m' (vh) (r - Vh) + O {r - r^f (5.16) 

5 (r) = + 5' {tu) (r - r/,) + O (r - vuf (5.17) 



where / (r/j) and h (r/j) are now the shooting parameters, r/j is a free parameter, and the field 
equations (5.3)-(5.7) give the derivatives /' (r/j) , h' (r/j) , m' (r/j), and 5' (r/j) as functions of 
these parameters. 
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On the basis of past investigations [7-10], we expect solutions (for a particular choice of 
free parameters) to exist only for discrete values of the shooting parameters. Using the above 
analysis, we can anticipate the asymptotic behavior of the fields for a small neighborhood 
surrounding those discrete values in shooting parameter space. Since / exhibits oscillatory 
behavior for < 0, and > leads to diverging /, we should look for a parameter range in 
which I/I approaches unity at large r and either exhibits a turning point or exceeds |/| = 1 
and diverges. If for such a parameter range approaches and (for //^ nonzero) then under- 
goes a turning point and diverges, or becomes negative and diverges, then our neighborhood 
should contain a point which gives the correct asymptotic behavior. We explain the shooting 
procedure in more detail below. 

Although the vacuum values \f (oo) | — 1 and (oo) — are shared by all solutions, the 
behavior of the field equations as r ^ oo provides interesting distinctions between massless 
and massive dilatons. For massive dilatons, the leading-order asymptotic expansions of the 
fields and metric functions are 



f{r 
h{r 

m (r 
lnr(r 

S{r 



± 



(-;) 



~ ae 



-m4,r 



M 
In 



1 



+ 



M 



1 c^e-^-^"" 

-(^0 + 77- 



(5.18) 
(5.19) 

(5.20) 
(5.21) 

(5.22) 



2 r-^ 

where c and a are positive constants, Tq and 5q are the rescaled metric constants introduced 
in the discussion of (2.18), and the +(— ) sign in (5.18) corresponds to an even(odd) number 
of nodes in the function / (r). Note that the presence of m^, defined by 
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(5.23) 



forces the dilaton to approach its asymptotic value exponentially as we expect. For a^o ^ 1, 
we find it convenient to use in place of //^ as a free parameter. In the massless case 
jj? — 0, the expansions are 

2/ioo - c) " 



/(r)~± 



c 1 c 6M 
-1 + - + -^ 



(5.24) 
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/i(r)~/ioo + ^— + g -2 (5.25) 

m{r)^M--^- (5.26) 
In T(r)~ In + (5.27) 

where we have included the first two nontrivial orders in 1/r to help demonstrate that the 
constant hoo plays the role of a Coulombic charge [49] in one metric function 

RHr)^{^-?^.'^.o(l)y\ (5.29) 

which has the form of the Reissner-Nordstrom solution, but is absent from the unrescaled 
form of the other: 

l-^ + ffi)V (5.30) 



2^2 y J, yj.. 

Noting that hoo — hm ( —e~^r^<^' I i?'^ ) , we can obtain an integral expression for h^o [49] 
from (5.2), 



ro HI \ / Jro 



2p-2{h/r+<t>o.) ( lL_^^ _ r-V 



(5.31) 

We have included V' in this expression for the sake of generality, and have denoted the lower 
limit ro to emphasize that (5.31) applies to both regular and black hole solutions, since 
(— e~^r20'/i?2) vanishes at ro = r/j as well as ro = 0. It is difficult to extract any further 
information about hoo from this expression alone, but we can derive a simpler expression 
for hoo and other useful relations by applying the scaling argument techniques of [53] and 
utilizing the metric relation (3.17). Before doing so in the next subsection, we address an 
ambiguity which arises from the definition of m^. 

In the above discussion, we are careful to qualify "massless dilaton" by the condition 
jj? — V — because (5.23) indicates that there are two ways to achieve massless dilatons 
in this theory: m? vanishes if either //^ = or aoo — 1. The latter case, corresponding to 
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0CX) = 5^n(3/6o) ~ 0.830, is also special in another respect. If we consider a power series 
expansion of V about a = aoo, we find that the coefficient of the nth order term varies as 
(1/q;oo)^"'~^\ so ttoo = 1 niarks the radius of convergence of the potential expansion. This 
would present a restriction if we chose to work with the representation Y = e"^*^ for the 
dilaton preferred by some authors, but for our field choice the power series expansion is 
well-behaved for all cKoo > 0: 



2 

ttoo + 1 / (ttc 
^ 4 2 (4[l-llaoo]+«L["oo-3][/Q;oo->i>lJJ24^^o/'25\ 



f2 \ 4 2 ^«oo - 72 , 4 3 («oo " 1) ("oo - 3) 73 

V / V ttoo + 1/ (ttoo + 1) 

(4 [1 - llaoo] + «L ["00 - 3] [7aoo - 33]) 



(5.32) 



12 (aoo + 1)^ 
From this expression, we find 

V {cl>, soo) = |/ (0 - <Poof + ((0 - 0oo)') (5.33) 

for aoo = 1 and near = ^ln(3/6o), so the potential resembles that of a Xcf)^ theory 
with A = 3!/x^ in the vicinity of its minimum, while for all other positive values of ttoo the 
potential has (f) and (f) contributions. When we examine the asymptotic behavior of the 
field equations for this special case, we find the same leading order behavior as in (5.18)- 
(5.22) with the exception of the expansion for h (r), 

, , , 1 c2e~2<^°° 1 c^e-^-^- (2M + c) , , 

Hr)--^-^, 5 ^ ^'-'^^ 

Note that the behavior of the field is completely determined by the gauge field to leading 
order: although the presence of a cubic V' in (5.7) eliminates the possibility of nonzero hoo 
which arises in the /i^ = case, it does not otherwise influence the asymptotic behavior of 
h (r) until several orders beyond those shown in (5.34). Most significantly, in the ctoo = 1 case 
the dilaton field can only match the boundary conditions if < near r — oo, which from 
our discussion of (5.9) is inconsistent with the behavior dictated by the full dilaton equation 
for 00 > or (r/i) > 0. The same arguments also apply to the range (pmax ^ <pQ, <P ('"ft) < 0, 
where (pmax is the location of the local maximum of V (see fig .1). For the range (p < (pmax, 
where V' > 0, the dilaton field is not strictly forbidden to increase, but the exponential 
amphfication of the gauge field term relative to V' in (5.7) and (5.9) makes a solution with 
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increasing toward = highly unhkely. After checking this possibihty for a wide range 
of initial conditions, we believe that no solutions to EYMD theory exist for ctoo = 1- 

Further Analysis 

To further explore the EYMD system analytically, we make use of a scaling technique 
developed in [53] for regular solutions and extend it to study black hole solutions. We find 
that the same results may be obtained from the metric relation (3.17), but we investigate 
the equivalence of the two approaches elsewhere [71] . 

The general procedure of [53] involves defining a nonlocal energy functional 

M = j dri^e^^mj = j dre^'^C^^^ (5.35) 

where 

is a sum of (non-negative) contributions from the kinetic terms of the matter fields. 
Demanding that M be stationary with respect to variations of the independent matter fields 

then yields the Euler- Lagrange equations when we treat C^^^ and 6 as functionals of those 
fields, and we restrict our attention to theories for which the action may be written 

S(x j dtdr + e-\ (5.37) 

where U and V are independent of the metric functions. The combination C^^^ = m' — 5'm in 
(5.35) is just m' with the metric function R (r) set to unity; with natural units restored, 
is equivalent to an effective matter Lagrangian with G — 0. By introducing radial scaling 
transformations for the independent matter fields ip of the type 

(r) = A"*^/; (Ar) , (5.38) 

which lead to the decompositions 
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we can find the A-dependent energy functional M (A) and obtain tlie constraint 



dM 
"dX 



(5.40) 



Depending on the positivity properties of C^^^ and the various scahng constants the 
constraint can be used to cstabhsh non-existence theorems or useful virial-type relationships 
for the kinetic and potential contributions of the matter fields. In the present case, 



2r2 

V^ = r'V{ct>,s^), 



(5.41) 



and the constraint dM / d\\\=i = assumes the form 

) pOO pCO pOO 

dr {y'f + y'^)=2 dry} - 2 / dr 6fy' - 2 / dre 
Jo Jo Jo 

where Sf is the contribution to 5 from Uf, y = e~^m, and 



f 

Jo 



(5.42) 



y'f ^ e"' (Uf + Vf) y'^ ^ e-' {u^ + v^) 



(5.43) 



are the gauge field and dilaton contributions to y'. 

The approach to black holes in [53] involves defining an effective Lagrangian on 
a fixed black hole background and a different scahng for the matter fields, ijjxir/rh) = 
'ip{{r/rfi)'^). Such an approach yields a complicated integral relation which generically in- 
volves ln(r/r/j) in the integrand. Though still quite useful in establishing non-existence 
theorems, the result of the procedure is not very useful as a pseudovirial relation. Since we 
are interested in a means of simplifying integral relations such as (5.31), a relation analogous 
to (5.40) for black holes is more appropriate for our purposes. To extend the regular solution 
analysis, we write the energy functional as a sum of horizon and non-horizon contributions, 

/•OO 

M= / dre-'^£(°) + e-'^(''°)m (ro) . (5.44) 

Jro 

By treating , 5 and S (ro) as functional of the matter fields, we again recover the correct 
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equations of motion from the variation of M. The corresponding constraint equation is 



dM' 



A=l Jro 



g-5 _ e-^iro)^' (^q) 



(5.45) 



In the notation introduced above, this constraint can be written for EYMD theory as 
/ dr {y'f + y'^) =2 / dry'f -2 dr 5fy' - 2 / dre-% 

Jrn Jrn Jrn Jrn 



(5.46) 



2e-^(-o)^ (ro) 5f (ro) - roe-^(^°) {Vj (ro) + V4, (ro)) , 



the left-hand side of which is equivalent to M — e ^^^°)m (ro). 

With the addition of the definition x = — r^0', the integrated dilaton equation (5.31) in 
this notation becomes 



h 



°° d 
dr- 



ro 
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e I a; 

r 



2 dr [y'f - 5fy'] - / drr'^e-^V' + 2 / dr (Sfy)' . 

J ro J ro Jq 

(5.47) 

Substituting the constraint (5.46) then gives 



PCO POO 

hoo= dr {y'f + y'^) - dr r\-' {V -2V)+2 [6fy] ^ 

Jro Jro 

+ roe-^^'"^ {Vf (ro) + (ro)) + 2y (ro) 5f (ro) 

pOO 

=M - e-'^('^°)m (ro) - dr r'^e'^ [V' - 2V) 

Jro 

+ roe-'^(''°) (y^(ro) + T/^(ro)). 



(5.48) 



With jj? — 0, this establishes the simple result /iqo — M for regular solutions and 



hoo^ M + rue 



(5.49) 



for black holes with ro = r/j = 2m (r/j). 
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For nonvanishing potential, eq.(5.48) also gives the interesting regular solution relation 



M 



(5.50) 



where e^'^Vgtr = V and "str" describes quantities in the string frame, which is related to our 
Einstein frame by the conformal transformation g^^ = e~^'^g^^ : 



J d^x./^re^'^^ {Rstr + 49^09^0 - Wstr (0, s) + . . .) 
I — > J d^xv^ {R - 2df,(f)d^'(l) - W (0, s) + ...). 



(5.51) 



If there is any significance to the suggestive form of (5.50), it is not readily apparent. 

The regular solution scaling analysis of [53] and its extension to black hole solutions 
provide an analytical relation among the solution parameters and the integrals of the dilaton 
potential which can be verified numerically. Since this nontrivial information about the 
system must be realized by the field equations, we expect that it might be obtained directly 
from them. A closer examination of the metric relation (3.17) reveals this to be the case. In 
the notation of the scaling discussion, (3.17) may be rewritten 
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2y'f + 2yS'f-2e-%. 



(5.52) 



Using the T equation 



i?2 



m 
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(5.53) 



to express the boundary terms, we obtain the relation 



M + 



(InT)' 



r=ro 



--M - e-'^('-°)m (ro) + roe-'^(^°) {Vf (ro) + (ro)) 

pco poc r-oo 

--2 / dry'f + 2 / dr S'^y - 2 dre~^V^, 

J To Jrn Jrn 



(5.54) 



which is identical to the scaling result (5.46) with the 6fy' contribution integrated by parts. 
This result provides the alternate expressions 



hoc + 



(InT)' 



(5.55) 



r=rh 
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M 



9str e 



-2(j> 



str 



(InT)' 



(5.56) 



r=rh 



for eqn.(5.49) and eqn.(5.50) generalized to black hole solutions. Note that the different 
signs of the kinetic contributions in m' and r (^uT)' / F? indicates that the result from the 
functional approach (5.45) agrees with (5.54) only when tq is the radius of an event horizon, 
not some arbitrary radius. It is also interesting to note that the regular solution parameter 
values correspond to the limit of the black hole values, which is consistent with solu- 

tion properties observed in previous studies of non-Abelian gauge fields coupled to Einstein 
gravity [7-9]. 

We can use the combination of the metric relation (3.17) and the dilaton field equation 
to determine an additional analytical relation. The combination 



d 
dr 



1^(0' -(InT)') ^r\-'{V'-2V) 



(5.57) 



can be integrated with the help of (5.48) to yield a useful relation between (j) and T: 



(0(r)-</'oo)=lnT(r) + (M-/ioo) / d?- 

Joo 
J oo 



(5.58) 



+ / dr- 



In the iJ? — case, we have the simple regular solution feature (r) 



In T (r) , which 



can be easily verified numerically, while for black holes the above equation gives 



(0 (r) - 0oo) = In T (r) + (M - /loo) / dr- 

J oo 



(5.59) 



The generic relation for a nontrivial potential is obtained from (5.58) by setting /iqo to zero. 

Numerical Regular Solutions 

As we observed above, solving the the EYMD equations numerically is a two-parameter 
shooting problem. In terms of the shooting parameters b = f" (0) /2 and 0o = ^ (0) = 
{h' (0) -I- 0oo) the boundary conditions (5.10)-(5.13) become 



/ (r) = - 1 + fer^ + —b (-24b^e-^'^° - 96 + 12Vq + AV^) ^ O (r^) 
30 V ^ 



(5.60) 
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InT (r) 



2m (r) 



/i(r) 




(5.62) 



(5.63) 



(5.61) 



- —626-2-^0 C28862e-2<^o - 1446 + 144Vb + 48^0') + C» (r^) 



where Vq, Vq' ^'^'^ ^o" potential and its derivatives with respect to at = 0o, and 

the shooting parameters satisfy b > and 0o > 0oo- We evaluate the initial conditions 
a,t r — 10~^ and use global error tolerance 10~^2 an adaptive fifth-order Runge-Kutta 



toward r — oo. For a range of b and 0o in the vicinity of a solution, the fields behave much 
as we anticipated: the gauge field either undergoes a turning point at |/| ^ 1 or diverges, 
and (f) either undergoes a turning point at ^ or becomes negative and diverges. In 

the /^2 = case, the property (p' < excludes turning points, and we instead observe 
a transition between monotonic decrease to (cxo) = and the divergence —oo over 
a small parameter range. In terms of h, the finite dilaton behavior is characterized by 
a turning point after exponential decrease toward h (oo) = 0, or the approach toward a 
positive constant value hoc when fj? — 0. 

In the massive case, finding a neighborhood of some discrete point (6, 0o) which exhibits 
these properties does not guarantee that one has found a legitimate solution to the EYMD 
system. The key to determining whether the results of the shooting procedure constitute 
valid solutions lies in the exponential behavior of h. For h to decay exponentially to zero, the 
gauge field coupling term in (5.7) must be insignificant. Since this contribution is positive 
definite at finite radius, we must rely on its algebraic approach to zero via (5.18) to satisfy this 
condition. As fj? is increased from zero, the radius at which h exponentially decays behaves 
roughly as r ~ l/m^, and eventually encroaches upon the fixed region where the gauge 
contribution algebraically decays. In other words, the screening of the Coulombic dilaton 
charge hoo occurs at a radius r ~ l/m^ which approaches the zone where the local magnetic 



ordinary differential equation solver, adjusting (6, 0o) for fixed 0oo a-nd and integrating 
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charge density vanishes. According to (5.7), once is large enough that these regions 
overlap, the gauge field source drives the dilaton away from its vacuum value and solutions 
are not possible. For a range of near this overlap, h appears to decay exponentially 
and satisfy the numerical solution criteria, but close examination reveals deviations from 
exponential behavior which prohibit extrapolation to h{oo). To determine a maximum 
allowable mass {mfj))max in practice, we must set a limit on the deviation of h from the 
behavior required to match the boundary conditions. By considering next-to-leading order 
terms in (5.7), we find that 

h - ar-^'^^ e-^'^' (5.64) 
describes the asymptotic behavior of the dilaton more precisely than (5.19), so that 



5h" _ 1 



h"-ml{l + —\h 



(5.65) 



gives a fair measure of the deviation from ideal behavior. Another useful quantity is 



Ah" 
h" 



h" 



,-2(Vr+0oc) 



h 



2r2 



(5.66) 



which directly measures the contribution of the gauge field coupling term to h" . In the asymp- 
totic regime of a valid solution, we expect \A.h" /h"\ -C /i" <S 1 and (5.66) to be comparable to 
(5.65), but the maximum acceptable value of either quantity is somewhat ambiguous. Since 
the size of the contributions to the dilaton equation which are not accounted for by the de- 
viation formula (5.65) are roughly of order h" /r^, and f ~ 10^ is the characteristic radius at 
which / obeys the asymptotic form (5.18), we adopt the criterion \Sh"/h"\ ^ 1/f^ 10~^. 
We find that this criterion gives the consistent result (m0)max 
of solutions we investigate. 



1/r ~ 10 ^ for the classes 



When actually obtaining solutions, we adjust the the shooting parameters until the 
solution bracketing conditions (the turning points and divergent behavior of / and h which 
characterize the solution neighborhood) indicate that the intervals containing the discrete 
solution values are smaller than our machine accuracy. To achieve this precision as is 
increased, we truncate h by taking /i — > or by attaching an exponential tail h — > ae"'"'*^ 
at the final turning point; this allows the integration to proceed so that the / bracketing 
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condition can be determined. To justify this procedure, we apply the deviation criterion 
(5.65) at the turning point and verify that h behaves according to (5.64) to better than one 
part in 10^. Since h tends algebraically to hoo rather than experiencing a turning point when 
the dilaton is massless, the truncation procedure is unnecessary and the numerical pitfalls 
posed by the final term in (5.7) disappear. 

The results of the shooting procedure for the choice 0oc = are displayed in fig. 2, with 
the massless and {m^)rnax solution properties summarized in table 1. Like the results of 
previous studies [7-10], solutions can be classified by the number of nodes k exhibited by 
the non-Abelian gauge field function /. Though an infinite number of solution classes exist, 
we focus our attention on the lowest odd- and even-A; classes. For both of these classes, 
we performed an identical shooting procedure with the toy potential V {(f)) — m^(l)^/2 for 
comparison. As table 1 indicates, the results agree to better than one part in 10^ for the 
narrow range of allowable dilaton masses, which indicates that the higher-order terms 
in the potential expansion (5.32) are negligible for this choice of 0oc- Though the total 
mass M measurably increases as k increases and varies over the allowed range, the only 
substantial change in the function plots occurs for h (r) . It exhibits a maximum at small 
radius and approaches hoo at large r, where the Coulombic dilaton charge is exponentially 
screened progressively closer to the decay zone of the gauge field. The plot of the actual 
dilaton field h (r) /r, along with the gauge function / (r) and the mass-energy m (r), exhibits 
nontrivial variation only in the decades surrounding r — 1: the characteristic radius of the 
solution is fixed by the string coupling, which in our dimensionless variables mimics the 
choice 'g = 2/\/a' = 1. Note that the maximum dilaton mass, which we only determine to 
one decimal place using our imprecise criterion, decreases as the characteristic radius of the 
gauge field decay zone increases with increasing k. 
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3 X 10- 
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8.3620815 


1.7927935 
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0.6848332 


1.7927935 
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8.3391611 
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1.7928245 


2* 
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-3 


8.3395171 


1.7914400 


0.6853065 





1.7928145 



* (run with V ((f)) — ^i^cf? for comparison) 
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Table 1. 



The solution parameters M, hoo and InT(O) in table 1 are determined with the aid of 
the asymptotic expansions (5.18)-(5.21) and (5.24)-(5.27). In the massless case, they provide 
a good check of some of the analytical relations determined above. In particular, hoc — M 
to better than seven figures in accordance with (5.49), while the prediction of (5.59) that 
(f) (r) = In T (r) is confirmed at r = (and other points) to an accuracy exceeding the global 
error tolerance 10~^^. The agreement of numerical and analytical results is compelling 
evidence for the general accuracy of our shooting method; the numerical results of other 
authors who recently considered the massless case [49] do not clearly exhibit these relations, 
though they are in general agreement with table 1. 

As values of 0oo in the range —0.8 ^ 0oo ^ 0.8 are used in our procedure, solutions 
appear to be related to the 0oo = solutions by the scaling of the radius and some physical 
parameters. To understand and quantify the scaling, we consider the simplified case of the 
massless dilaton. When we ignore the axion and the dilaton-axion potential in the action 
(2.14), the dilaton explicitly appears only in the exponential coupling to the non-Abelian 
field strength. If we absorb the constant 6"^*^°° into the gauge coupling and then rewrite 
the theory in terms of dimensionless variables and parameters, we recover the 0oc = 
theory but with 'g e~'^°°'g. From the definitions of the dimensionless quantities (3.2)- 
(3.3), we therefore expect the radial structure of the 0oo solutions to scale according 
to r — > e~'^°°r. The mass-energy should similarly scale as m — > e~'^°°m, but the amphtude 
of the dilaton deviation field and the gauge field should remain unchanged. It follows 
that the dilaton shooting parameter shifts such that (0) = (pQ — (poo is unaffected, while 
the gauge field parameter scales as h ^ he^'^°° to compensate for the scaling of the intital 
value of r^. Though the introduction of the potential and its complicated dependence on 0oo 
changes this picture, solutions approximately obey the same scaling relations for the range 
of 0OO examined, with ( 

n^(j))max — ^ ^ {j^(j))max ^ expected from (3.2) and the relationship 
{nT'(j))max ~ 1/f discussed above. As 0oo approaches the critical value ^ ln(3/6o) ~ 0.830, at 
which the potential (5.32) reduces to the form (5.33) to leading order, solutions become 
harder to obtain numerically and appear to be forbidden at | ln(3/6o) for reasons examined 
above. For 0oo above this critical value, the leading order term is restored in the 

potential expansion and solutions are again possible, though we do not explore the scaling 
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properties of this solution region in depth. 
Numerical Black Hole Solutions 

To find numerical black hole solutions, we use the conditions 



/ (rh) 
h' {rh) 
m! (rh) 

(Th) 



[1 - f {rh)] f {rh) 



rh - e-^^('>^)[l - p {rh)f /rh - 2V (0 {rh)) r^) 

V (0 {rh)) rj - e-^^(^^) [l - f (r,)] ^ /r, ^ h (r,) 
(rh - e-20(-.)[l - /2 {rh)f /rh - 2V (0 {rh)) r^) 



,-20m [1 fHrh)] +v{<l>{rh))rl 



h' {rh) 



h{rh) 



rh 



(5.67) 

(5.68) 

(5.69) 
(5.70) 



on the horizon, and use / (r/j) and h {rh) = rh (0 {rh) — (poo) as shooting parameters for 
rh — 1- The asymptotic properties of the fields, which we use to locate the neighborhood 
of a solution in shooting parameter space, are identical to the regular solution properties, 
so we follow precisely the same shooting procedure detailed above. In particular, the same 
truncation of h and maximum dilaton mass criterion are used to determine solutions in the 
massive case. 
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-.1321129 


0.5444857 


0.8651225 





0.1510935 


2* 


3 X 10-4 


-.1321128 


0.5444860 


0.8651223 





0.1510932 



(run with V {4>) 



)^ for comparison) 

Table 2. 



We again examine only the k — 1- and k — 2-node solution classes for the choice ^oo = 
and compare results for the full potential (2.15) and the toy potential V {(f)) — m^0^/2; the 
results are shown in fig. 3 and table 2. Again only the function h varies significantly as 
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TTifj, increases over the small allowed interval, with h approaching hoc and then vanishing 
exponentially at progressively smaller radius. The solution structure is nearly identical to 
the r > 1 portion of the regular solutions, which in part reflects the occurrence of the event 
horizon at the characteristic radius r = 1 of the dimensionless system. A closer examination 
of the regular solution functions reveals a sharp peak in the metric function R (r) near r — 1, 
where 2m (r) / r closely approaches unity, thus demonstrating that even the regular solutions 
are strongly gravitating. As we might expect from previous work with black hole solutions 
to theories with non-Abelian gauge fields [7-9] , the black hole solutions reduce to the regular 
solutions for the same choice of {(pooi 

rriff,} in the hmit r/j ^ 0. 

Though the total mass-energy and dilaton charge do not obey the simple relationship 
hoc — M enjoyed by the massless dilaton regular solutions, the analytical result (5.49) does 
relate /iqo to M, Sq — 5 (r/j) and the integration parameters r/j, /(r/j) and (f){rfj). Again, 
the results of our numerical procedure verify a nontrivial relationship between the physical 
parameters of solutions to better than seven significant figures. Black hole solutions for the 
massive dilaton also exhibit the approximate scaling properties explored above for non-zero 
0OO) but the specific relations for r, m (r) and {mfp)max only hold when we scale the horizon 
radius according to — > e~'^°°r/j. 

6. Conclusions 

In this paper we have studied static, spherically symmetric regular and black hole so- 
lutions to SU{2) gauge theory coupled to a massive dilaton, massive axion, and Einstein 
gravity. Our intentions have been two-fold: to explore solutions in the physically relevant 
context of low-energy string theory with massive scalar fields, and to determine whether 
"stringy" scalar fields lead to non-Abelian solutions with primary hair and good prospects 
for stability. After analyzing all the possibilities for fundamentally non-Abelian solutions, 
we found strong numerical evidence for regular and black hole solutions of a massive dilaton 
coupled to the Yang-Mills field (EYMD theory), and established a deeper understanding of 
certain solution existence techniques [53] in the course of exploring the solutions analytically. 
Though the case of a massive axion coupled to the gauge field appeared promising, we found 
that the full theory, which describes a massive dilaton and massive axion coupled to a dyonic 
non-Abelian configuration, is the only other situation which can admit solutions. We pre- 
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sented no numerical evidence for such solutions, but we were able to construct a consistent 
solution scenario. 

An important issue that we have not addressed in depth is the stability of our solutions. 
As we noted above, the primary hair solutions to EYMD theory are structurally very similar 
to the solutions of EYM [10, 7-9] and EYMH [23] theories, which have been interpreted 
as generahzed sphalerons and are generically unstable. Though examples of stable solu- 
tions with non-Abehan structure have been found, including "black holes inside magnetic 
monopoles" [20-22] and Skyrmion black holes [27-28], such solutions typically possess a net 
gauge or topological charge (and the resultant imprint at spatial infinity) which is lacking in 
our solutions. These observations, and the linear analysis of [46] which established the in- 
stability of EYMD solutions for a massless dilaton, make the stability of our massive dilaton 
solutions very unlikely. 

In light of this conclusion, one might question the relevance of pursuing numerical so- 
lutions to the full theory. As the only examples of gravitating SU (2) solutions with both 
magnetic and electric fields, such solutions would be interesting in their own right, but the 
lack of a net electric charge (which follows here from the asymptotic behavior of the field 
equations) would appear not to improve the chances of stability. It is conceivable, however, 
that the structure arising from the coupled electric and magnetic charge densities substan- 
tially modifies the sphaleron character of EYMD solutions, even in the absence of a net 
charge. Only in such circumstances, it seems, could we reasonably hope for stable solutions. 
Since it requires a four-parameter shooting procedure, the task of finding such solutions 
could present enough obstacles that these questions might remain unanswered. Though 
we have not yet attempted to obtain such solutions, we have described a strategy which 
might simplify this formidable task. We hope to test the efficacy of this strategy in future 
investigations of string-inspired non-Abelian dyonic solutions. 
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FIGURE CAPTIONS 

1) The dilaton-axion potential (2.15), which is of the form used in the study of SUSY- 
breaking via gaugino condensation in string theory. In (a), the rescaled axion field 
s = s(3/47r6o) is fixed at one of the degenerate minima s = Sqo = n for integer n, 
and the dilaton at spatial infinity is chosen to be (poo = 0. Solutions to Einstein- Yang- 
Mills-Dilaton theory correspond to rolling monotonically to the minimum from the 
right, confined to a region where V is well approximated by leading-order (0 — 0oo) 
behavior. In (b), the dilaton is fixed at 4> — 4>oo and the potential assumes the form 
V oc 2 (1 — cos(27rs)). A possible solution scenario for the full theory involves the 
axion traversing one of the maxima in the s direction and monotonically approaching 
an adjacent minimum as — > 0oo from above. The fact that such non-Abelian dyonic 
solutions are conceivable critically depends on the presence of both the massive axion 
and massive dilaton fields. 

2) One- and two-node regular solutions to Einstein-Yang- Mills- Dilaton theory for the 
dilaton potential of fig. la. The connection function (1 -|- /(r)), total mass-energy 
m (r) , dilaton (p (r) and h (r) = r0 (r) are plotted as functions of radius for a range 
of dilaton masses < m,^ < ("^0)^^^. The maximum mass ("^0)^^^ is reached as 
the exponential screening of h (r) , which asymptotically resembles a Coulombic charge 
hoo = h (00) in the m<^ = case, occurs near the zone where the local magnetic gauge 
charge vanishes. The other functions do not vary appreciably in this mass range. The 
exact correspondence between hoo and M = m(oo), as well as (r) — 0oo and the 
metric function In T (r) , is proven analytically in Section 5 for the massless dilaton 
case. 

3) One- and two-node black hole solutions to Einstein- Yang-Mills-Dilaton theory for hori- 
zon radius — 1 and the dilaton potential of fig. la. As in the regular case, only 
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the dilaton function h (r) = rcj) (r) varies appreciably as the dilaton mass approaches a 
maximum value. Though and M are not equal for = 0, an analytic expression is 
derived in Section 5 which relates these quantities to other global solution parameters. 
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